Topological quantum phase transition in the BEC-BCS crossover phenomena 
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A crossover between the Bose Einstein condensation (BEC) and BCS superconducting state is 
described topologically in the chiral symmetric fermion system with attractive interaction. Using 
a local Z2 Berry phase, we found a quantum phase transition between the BEC and BCS phases 
without accompanying the bulk gap closing. 
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Introduction- Recent progress in ultracold atomic 
Fermi gases, the Feshbach resonance which controls 
strength and sign of effective interaction, realizes the 
Bose-Einstein Condensation (BEC)-BCS crossover—""— — 
Bose condensation of bosonic molecules which are real- 
ized in real space pairing of fcrmions and the BCS su- 
perfluid with Cooper pairing in momentum space. A 
long time ago, Leggett^ proposed a trial wavefunction 
for the ground state which describes not only weakly at- 
tractive case where formation of the Cooper pairs occurs, 
but also the strong attractive case where a dilute gas 
of molecules undergoes the BEC. Also substantial num- 
ber of theoretical works on the BEC-BCS crossover have 
been performed^""— after this pioneering work. By using a 
Ginzburg-Landau type argument based on the symmetry 
breaking and local order parameter, these states are not 
clearly distinguished. They are considered as a crossover, 
that is, the two ground states are adiabatically connected 
even in the thermodynamic limit. 

Both of the BCS and the BEC ground states are typ- 
ical examples of gapped quantum states, which we call 
quantum liquids, in the sense that there are no symme- 
try breaking phase transition between them. Still they 
are quite different and characteristic quantum states. As 
discussed in this paper, there are clearly distinguished by 
using a topological quantities. Therefore they belong to 
a new class of matter as the topological insulator based on 
a novel classification scheme (topological classification). 

A class of quantum Hall states is a typical ex- 
ample of topologically non trivial quantum liquids 
where topological quantities play fundamental roles for 
the characterizatiori^""— . Also a time reversal invari- 
ant analogue of the quantum Hall states are studied 
intensivel y ^^1^^ . 

As for the quantum liquids with non trivial topological 
structure, we are proposing to use topological quantities 
such as the Berry phases and the Chern numbers using 
Berry connectionsii""— . 

An important characteristic feature of the topologi- 
cal insulators is an appearance of local degrees of free- 
dom near boundaries and impurities as generic edge 
states^iiiii^ii^. Although the bulk is gapped and fea- 
tureless, the edge states characterize the topologically 



non trivial bulk. They are not independent but in- 
timately related each other, which is known as hulk- 
edge correspondence, where topologically non-trivial bulk 
guarantees the existence of localized modes and such 
low energy localized excitations characterize the gapped 
bulk insulator— conversely. The quantum Hall state of 
the graphene also belongs to this topological insulator 
where the bulk-edge correspondence is important for the 
description22r— . The bulk-edge correspondence is also 
realized as the existence of the Kennedy triplet for an 
open integer chainS^. Such characteristic edge modes ap- 
pear in the valence bond solid (VBS) states^! and spin 
ladder with cyclic exchange interactio n^^i^^ . 

In this paper, we introduce a local U{\) twist for the 
Hamiltonian and define a Berry phase using its many 
body ground state, which is quantized into Z2 due to the 
chiral symmetry of the Hamiltonia n^^'^^ . Although the 
bulk gap is adiabatically connected between the BEC- 
BCS crossover, the gap of the twisted system may col- 
lapse at some value between the BEC-BCS crossover. It 
is allowed since the gap is not a thermodynamical prop- 
erty and can be collapsed by the local perturbation. It 
actually occurs in the BEC-BCS crossover. Then the 
crossover of the bulk is distinguished by the local quan- 
tum phase transition. 

Model and Zi Berry phase - Let us start from the at- 
tractive Hubbard model at half-filling and discuss a mean 
field Hamiltonian 
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where q.o- is annihilation operator at site i with spin a. 
The number of sites is N. The summation over < i,j > 
is restricted to the nearest neighbor pairs. The density 
with CT-spin at site j is given by Uj^^ = cj ^cj^o-. We im- 
pose the periodic boundary condition. We assume the 
lattice is bipartite in 1 to 3 dimensions. To make the 
discussion clear, let us perform the following particle- 
hole transformation in the mean field hamiltonian Hmf, 
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FIG. 1: Lattice structure of the efTective mean field system 
for the two dimensional case. We introduce a twist on the 
bond defined in Eq.((2]). 



The order parameter of the superconductivity A — 
— \U\{cl^cl ^) = —\U\{uldi), is chosen as a non-negative 
real number. In this paper, we do not perform the 
self consistent calculation and just assume that the or- 
der parameter A is finite. Here the expectation val- 
ues (O) is taken for the ground state. The one-body 
eigenvalue problem of the iV-site system is given by 
-HuFlq^^'^) = -B(J) < < ... < £;(2W)), 

where the one-particle eigenstate is l^^-'^) = ((/)^-'^)^|0) = 



) with an orthogonal nor- 
malization condition q '^m"''*0m"'' ~ Sij- Here the 
vacuum state |0) is defined as Uj|0) = = dj\0) for 
any site j. The Af -particle eigenstate is constructed as 

|0)m = n,^i(0^'^)^|O) with the eigenvalue EJIi^^^'^- 
At half-filling band, by Fourier transform it is easily 
found that the ground state \(I))n has a finite excitation 
gap energy 2 A for finite interaction U for any bipartite 
system. 

The BCS and BEC states are not clearly distinguished 
and they are adiabatically connected as a crossover even 
in the thermodynamic limit. In order to define the Berry 
phase, we modify the terms locally only at j = a in Eq.([T|) 
as follows, 
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Hereafter we take t as unity. Figure [2] presents the 9- 
dependence of the energy spectrum of Eq.([3]) for the 
square lattice. The energy gap 

minimizes at = tt. In Figl2] (b), we can see the level- 
cross occurs at {0,E) = {"^,0) around A = A^ ~ 1.250 
for two dimensional case. As we will see later, A = Ac 
is the critical point between the BEC and BCS phase of 
the BEC-BCS crossover, although the gap of the trans- 
lational invariant system {9 = 0) is always open. 

Due to this local perturbation, the spectrum changes, 
which can be regarded as the edge mode. It characterizes 
the gapped bulk feature as a realization of the bulk-edge 
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FIG. 2: Energy spectrum of Hmf(^) for the two dimensional 
case {N = 16x 16): (a) A = 0.5. BCS phase, (b) A ~ 1.250. 
The level cross occurs at {9,E) = (7r,0). (c) A = 2. BEC 
phase. 



correspondence. Under the above deformation, the half- 
filled eigenstate \4))n depends on the 6, which denotes 
|0(6'))jv. The Berry phase 7 is then defined as 
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deN{m\9em)N- 



(2) 



The transformation {uo,d,) {—d,,Uo) (see Fig. [T]) for 
one sublattice reduces the Hamiltonian in the following 
bilinear form: 



(3) 



This Hamiltonian ^ has a chiral symmetry, 
{Hmf,^j} = 0, 7^ = 1. Unless the determinant 
of the matrix Dn{9) vanishes, the half- filling ground 
state of the Hamiltonian (jSj has a finite gap. Therefore, 
the Berry phase 7 is quantized as or tt modulo 27r (Z2 
Berry phase)^. Using this Z2 Berry phase, we identify 
whether the half-filled ground state is a BEC or BCS 
state. The Berry phase ^ for the half filled ground 
state of the chiral symmetric Hamiltonian can be 
obtained as following exprcssiork^: 



7 = 



d9 ImlogdetDAr(6'). 



(4) 



By definition, we can show that the Z2 Berry phase re- 
mains invariant through an adiabatic deformation until 
the level cross occurs. That is, the Berry phase is topo- 
logically protected. In this paper we evaluate the Z2 
Berry phase intuitively using an adiabatic continuation 
without direct calculation of Eq. (|4]) . 

Adiabatic deformation- We consider the two types of 
the adiabatic deformation to understand the BCS and 
BEC states explicitly. Since we introduce the flux 9 on 
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(a) BCS regime 



(b) BEC regime 



FIG. 4: Two types of the adiabatic deformations. One defor- 
mation is the decrease of the local bond (type (i)) (a). The 
other is the decrease of the hopping connecting the deformed 
local bond being zero (type (ii)) (b). 



one local bond, it turns out that det Dn{9) has a form 
Ae'^ + B (A and B are real and independent of 6) by 
the Laplace expansion. For a finite attraction U (i.e., 
finite A), only at 6* = tt, detZ)jv(0) can become zero 
where the level-cross occurs at zero energy. During the 
change of the variable 9 from to 27r, the determinant 
detI?Ar(0) draws a closed curve in the complex plane as 
shown in FiglS] If the det Dn{9) winds around the origin 
of the complex plane m times when 9 varies from to 
27r, the Berry phase 7 is given by 7 = nm modulo 27r. 
When the coupling strength satisfies that < A < Ac, 
the closed curve of det Dm {9) does not enclose the origin 
(see FiglSl^a)). From Eq.(|3]), we obtain the Berry phase 
7 as 7 = 0. When A > Ac, the closed curve det£>Ar(0) 
winds origin once (see Fig|3](b)). Then, we obtain 7 = tt. 

Type (i) : As seen in Fig. ^a), for small A(= 0.5), 
det Dm (9) draws the circle without surrounding the ori- 
gin in the complex plane. Then we adiabatically decrease 
the magnitude of the twist interlayer hopping A at j = a 
to zero (see Fig. |4] (a)). During this adiabatic deforma- 
tion, the circles in the complex plane shrink to the one 
point in a concentric fashion without touching the ori- 
gin as the interlayer hopping A on the twist bond goes 



FIG. 5: Relation between the s-wave gap A and the attractive 
interaction strength \U\ in one dimension. In weak coupling 
(BCS) regime, A behaves A ~ 8 exp — 2/(7r|f/|), where the 
Berry phase is 0. In strong coupling (BEC) regime, A ~ 
\U\/2 — 2/\U\, where the Berry phase is vr. 



TABLE L Critical values Ac for dimensions d = 1, 2 and 3. 



Dimension d 


Ac 


d = 1 (linear chain) 
d = 2(square lattice) 
d = 3(cubic lattice) 


2/^~ 1.1547 
1.250 
1.588 



to zero (i.e., during the adiabatic deformation ^ ^ 
and B remains same when we write down deiDN{9) as 
det Dm{9) = Ae^^ + B). During this deformation, the ex- 
citation gap of the half-filling ground state remains finite. 
After the deformation, the dependence of 9 disappears in 
Eq.([3]). Due to the invariance of the Z2 Berry phase un- 
der the adiabatic deformation, the Berry phase 7 should 
be zero for small A, where the binding of the pairing 
is in momentum space (BCS state). In contrast, dur- 
ing the deformation, the circle touches the origin in the 
complex plane once for large A(= 2.0). For large A the 
level-cross occurs at zero energy for large A. At the level 
crossing point during the deformation, the Berry phase is 
not protected any more. To understand the Berry phase 
for large A we consider another type deformation. 

Type (ii): We decrease the hopping t connecting the 
site a into zero. Figure [HJb) shows the detDAr(6') draw- 
ing the circle with surrounding the origin once in the 
complex plane. In this deformation, the circle is parallel 
shifted so that the center of the circle becomes the origin 
of the complex plane i.e., A remains same and i? — > 
for detD]y{9) ~ Ae^^ + B. During this deformation the 
level-cross does not occur at zero energy for large A. Af- 
ter the deformation, the twist interlayer hopping A at 
j = a becomes isolated. Then, it is easy to show the 
isolated interlayer bond gives the tt Berry phase. Thus, 
the two electrons form a dimer bound state (BEC state) 
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for large A. When A is small, the Berry phase is zero as 
discussed above. During the deformation of type (ii) , the 
circles touch the origin in the origin, that is, the level- 
cross occurs at zero energy. From numerical calculation 
we find that there exists only one critical point A = Ac, 
where the Berry phase is (tt) for A < Ac (A > Ac). In 
a similar manner, the critical value Ac for the other di- 
mensions d can be calculated (see Table It is obtained 
numerically as Ac ~ 1.588 for the cubic lattice. 

In one dimensional case, the expression of detDjv(^) 
can be obtained analytically as 



deti:>w(6') = iAe'^UN^i ( ^ 



, iA 

2Un-2 ( ^ 



Here the J7„ [x) is the Chebyshev polynomial of the sec- 
ond kind^l, defined as C/„(cosa;) = sin[(n -I- l)x]/sinx. 
Parametrizing A by A = 2sinha, the asymptotic form 
of the dct Dff{9) for large N can be expressed as 



det Dn{0) 



iAe 



2 sinh a 



sinh( 



cosh ( 



We have a real root on a for the equation det Dn{9) — 0, 
i.e., e'^ sinh a + e^" = only when 6 = tt. The solution 
is obtained as a (log3)/2, i.e., Ac = 2/\/3 ~ 1.1547, 
although the bulk gap {9 = 0) remains finite. 

For one dimensional system, the open boundary chain 
of Hubbard model are analyzed22i2^, which might be rel- 
evant to characterize the bulk features. 



Summary and discussions- We show the quantum 
phase transition in the mean-field attractive Hubbard 
model at half-filling. The Z2 Berry phase distinguishes 
the BEC-BCS crossover as a local quantum phase tran- 
sition, that is, the phases are separated by closing of the 
energy gap under the local twist, although the gap of 
the translational invariant system is always open. In the 
weak attractive interaction (BCS) case, the Berry phase 
is in one hand. On the other hand, in the strong attrac- 
tion (BEG) case, the Berry phase is tt (see Fig. [5]). That 
is, it characterizes whether the paired electron is itinerant 
or localized. This comes from the bulk-edge correspon- 
dence in the BEC-BCS crossover. Physically the itiner- 
ant Cooper pairs in the BCS phase are not affected by the 
local U{1) twist although the spatially localized bosons 
formed by the real space is crucially affected by the twist 
which results in the non trivial tt Berry phase. In this 
sense, our topological characterization reflects whether 
the size of the pairing is macroscopic or of the order unity. 
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